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Introduction  f ^ 

/)  ^ 

The  form  of  the  hulls  of  many  undersea  trail  is  approximated  hy  a 

closed  surface  of  revolution.  Tliis  property,  togctlier  with  the  numerical 

solutions  available  for  elastic  shells  of  revolution,  has  motivated  the 

study  of  these  figures  in  connection  with  the  .acoustic  radiation  problem. 

2 

See  [1-3]  . Solutions  have  been  given  which  involve  various  sim.pjifying 
assumptions  regarding  the  nature  of  the  radiated  acoustic  field.  [4-6]. 

Ref.  [4]  contains  a summary  of  some  of  the  recent  work  in  tliis  .trea. 

In  the  following,  the  scaled  acoustic  lield  equations  are  expressed 
in  circular  cylindrical  coordinates  and  the  circumferential  coordinate  sup- 
pressed through  a Fourier  series  expansion.  The  region  external  to  shell’s 
generating  a curve  is  mapped  conformally  onto  the  region  external  to  the 
unit  circle.  This  is  tantemount  to  the  construction  of  an  orthogonal 
coordinate  system  which  has  the  shell's  generating  curve  as  one  of  Its  coor- 
dinate lines.  The  numerical  conformal  mapping  technique  utilized  is  applicable 
to  infinite  regions  bounded  internally.  The  region  external  to  the  unit 
circle  is  mapped  into  the  region  internal  to  a rectangle.  The  cqviation  of 
motion  of  the  acoustic  fluid  together  with  the  shell  fluid  hoendarv  equations 
in  the  transformed  coordinate  system  is  expnssid  in  finite  difference  form. 

The  normal  shell  displacement  is  introduced  through  a ipeetral  1 1 pi  or-i'nt - 
ation,  as  a series  of  t lie  in  vacuo  .sliell  niodis.  Integrals  orruiiing  in  lor.ing 
terms  of  the  expansion  which  contain  the  unknown  .surface  aroui  tic  prei.sure 
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Abstract 


A numerical  solution  for  the  transient  vibration  of  an  arbitrary 
shell  o5  revolution,  surrounded  by  an  acoustic  medium,  has  been  formulated 
in  terms  ci  the  conformally  mapped  acoustic  equations,  fluid  sheli  boundary 
conditions  and  spectral  representation  of  the  shell  displacements.  The 
technique  readily  includes  the  boundary  condition  at  infinity,  internal 
shell  structure  and  is  exact  to  within  those  approximations  implicit  in 
the  finite  difference  method. 
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where  T = ct/s,  P = p/p^c^  and  = fluid  density.  In  (9),  the  angular 

coordinate,  x^,  may  be  suppressed  through  a Fourier  expansion  of  P 

rt> 

P(x,  ,x.,  ,x^  ,T)  = ) (P  cos  nx„  + P sin  nx„)  (10) 

n=0 

For  notation.-jl  simplicity,  P^  will  be  taken,  in  the  following,  to  represent 

either  P or  P . For  the  exterior  problem  x,  > 1.  However,  for  a closed 
nc  ns  1 — 

shell  of  revolution  there  exist  points  on  the  shell's  surface  for  which 
= 0.  The  corresponding  singularities  in  (9)  may  be  eliminated  through 
the  change  of  dependent  variable 


P = F /X 2 (l,x_) 
n n i J 


(11) 


The  region  external  to  the  unit  circle  = 1 may  be  mapped  into  the  interior 


of  a rectangle  through  the  inversion 


y = i/x, 


(12) 


See  [],c,3].  Substituting  (10),  (11)  and  (12)  into  (9)  and  taking  the 
Laplace  transform  with  respect  to  T,  assuming  zero  initial  conditions  gives 


3 ^ Fn  . , 3Fn  ^ _ , 3Fn  _ , 3'^Fn 


av 


‘1  3y 
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+ a:  — + a:  ^ + a!.  + a;  F^  = 0 (13) 
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'3  3 

where  = 1/y  - (DX^^/ 3x^)  / (y^X^)  , a^  = (3X^/3x^) /(X^y")  - 4 (3X^  (1  ,x^) /3x^)  / 

(y^Xj^d.x^) ) ; = 1/y^,  = (6(3X^(l,X2)/3x^)^/Xj^(l,X2)  - 2 ( 3^X^^  (1 , x^)  / 3^x^“' )) 

/(y^X^d.x^))  - 2(3X^/3x2)  (3Xj^(l,X2)/3x^)/(y^X^Xj^(l,X2))  - (g^^/y^)  (s^/a-^  + n‘/Xj  ) 
and  s is  the  Laplace  transform  variable. 

An  cfiei;Live  method  for  computation  of  the  coellicionts  occurriiig  in  the 
mapping  1 unction  (3)  is  given  in  [10].  This  technique  has  been  applied  to 
the  solution  of  the  cylindrical  problem  in  [7]. 

Eq.  (13)  m.ay  be  expressed  in  finite  diflerence  form  by  replacing  the 


spatial  deriv.it  ives  with  the  following,  diffctcncc  expressions: 


F, 


X x=x. 


Cn(x2>Fn,  n=l,3,  - 1/B,  = (a^-D/B  , = -u^/B.  B ’ •‘(a+l)h. 


h = X2  - Xi  and  ah  = x^  - X2*  For  points  not  adjacent  to  y = 0,  y = 1, 
x^  = -+-n/2  the  second  derivative  was  approxinated  by 

n=0,4  where  the  mesh  points  are  the  set  {x^,Xj^,X2,X2,x^} , ^ x^,  m=0,4. 


0,4,m^n,  D =xx,m=C,4, 
n mm 


in  n,  B = (A’^  - D )/2,  C = l(x  -x  ) , in 
n n n n n m 

m ^ n,  and  A^(X2)  = 12  X2^  - 6A^  x^  + 2B^)/C^.  For  points  adjacent  to  y = 0, 
y = 1,  x-  = tn/2  the  second  derivative  was  approximated  by  F,  1 = B (x  )F 

• j XX I X x^  n X It 

n = 0,3  where  = E x^,  m = 0,3,  m n,  = E(x^^-x^),  m = 0,3,  m ^ n. 


B 


(x,)  = (6x,  - 2A  )/C  and  the  mesh  points  consist  of  the  set  {x  ,x  ,x^,x  }. 
n 1 1 n n o 1 ^ j 

At  a general  nodal  point  y = y^,  = x^j^  the  finite  differerce  form  of  (13) 


is  then 


+2 


(D.^  , F.^  , + I>'.  , 

i+r,k  i+-r,k  i,k+r 


^,k.r^  - ° 


(1'.) 


r=-2 

where  the  coefficients  of  F in  (14)  are  found  by  substituting  the  above 
difference  expressions  into  (13)  and  the.  subscript  n has  been  omitted  in  the 
symbol  F^  for  the  sake  of  clarity. 


Boundary  Conditions 

At  the  shell  fluid  interface,  the  normal  fluid  velocity  roust  equal  the 
normal  shell  velocity.  This  implies  that 

aOP  /a  "n  + 31>  ,/3’n)  = -?-W  /3T-  (15) 

n nl  n 

where  P , P . and  W are  the  coefficients  of  the  Fourier  expansions,  in 
n tJ  n 

terms  of  X2»  of  the  scaled  reflected  and  radiated  acoustic  pressure,  the  scaled 
l.no'.Ti  incident  arrustic  pressure  and  the  scaleii  normal  shell  displacement 
respectively.  n is  an  exteriot  unit  vector  which  is  normal  to  the  shell's 
surface  in  the  unsealed  physical  coordinates  . Tt  nav  he  shown.  that 


9 P _ k ^ 
9 n ^ 9j^ 


(16) 


where  n are  the  rontravarient  components  of  the  unit  vector  n in  the 
coordinate  system.  The  tensor  transformation  law  implies  that 


k (ixb  -1 
n = n 

9x . 

1 


(17) 


Since  n is  normal  to  a surface  of  revolution  and  is  cylindrical 
_ 2 

it  follows  that  n =0.  This  and  the  fact  that  x^^  and  are  only  functions 
_ _ 

of  and  , togetlier  with  (17),  implies  that  n*"  = 0.  For  covarient  com- 
ponents (17)  becomes 


9xj 


- n . 

9x,  1 


(18) 


For  orthogon.il  coordinates  n =0  implies  that  = 0 and  therefore  n^  = 0. 
Summing  (18)  for  k = 3 gives  = (Sxj^/9x^)n^  + (9x2/9x^)n^.  However 

n = Vf/i7fj  = ((3f/9x^)ij^  + (3f/9x2)gj)/Vf  = n'*^  + and  n'  = , 

since  n^=  g^^n^,  = 0,  1 j,  g^^  = 1.  Here  ffx^^.x^)  = c 

is  the  equation  of  the  generating  curve.  Then  = ( (9Xj^/9x^)  (9f /9x^)  + 
Ox^/Sx^)  ( 9f /3x^) ) /jV£|  = (9f /9x^)  ^Vfl  = 0 since  fd.x^)  = c is  the  equation 
of  the  generating  curve  in  the  coordinate  system.  Tlien  n^  = 0 and  therefore 

3 

n =0.  Since  n is  a unit  vector  It  follows  from  the  definition  of  the 
scalar  product  and  the  definition  of  the  non-dimensional  coordinates  that 


Ik  2 2 3 1 2 

n n = l/a'.  Since  n “0  and  n = 0,  it  follows  that  (n  ) 

= l/(;!»gjj).  Substituting  into  (16)  gives 


1 or 


9 ri 


1__ 

■'"'hi 


(19) 


The  btnind.'iry  condition  (15)  Is  then 


J 


+ a 


(20) 


-1 


9F 


ilnL  = 

9 n 


_ ^ 
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where  y = l/x^^  and  is  replaced  by  through  (11).  The  finite  difference 
form  of  (20)  is  obtained  by  replacing  9Fj^/9y  with  the  first  difference 
approximation. 

In  the  rectangle  0 £ y 5.  1>  -n/2  i ’’^2,  (18)  must  be  satisfied 

for  y = 1.  The  boundary  at  infinity  corresponds  to  y = 0 on  which  “ 0* 

Assuming  that  Xj^(l,+  n/2)  = 0,  then  since  Pj^(y>  i ■”/2)  is  bounded,  it  follows 
from  (11)  that  F (y,  + t/2)  = 0.  Given  the  normal  shell  displacements,  V.'  , 
and  the  incident  acoustic  pressure,  the  boundary  conditions  together  with 

(10),  (11)  and  (13)  determine  the  reflected  and  radiated  acoustic  pressure 
over  the  region  external  to  the  cavity  defined  by  the  surface  of  revolution. 


Shell  Equations 

Solutions  of  the  shell's  equations  of  motion  may  be  introduced  into  the 
shell  fluid  interaction  problem  either  directly  through  finite  difference 
approximation  as  in  [2,3]  or  through  a spectral  representation  as  in  [6]. 

Since  there  exists  at  present  a more  comprehensive  library  of  shell  programs 
which  employ  spectral  representation  rather  than  integration  of  a finite 
different  approximation,  the  following  development  is  based  on  spectral  repre- 
sentation. For  notational  simplicity  U^,  U^,  and  W will  note  the  coefficients 
corresponding  to  some  given  value  of  the  summation  index  for  the  Fourier 
expansions,  with  respect  to  the  polar  angle  of  the  tangential,  axial,  and 
normal  shell  displacements  respectively.  Then  In  terms  of  generalised  cuorJi n.it e 

q (T)  and  scaled  normal  modes  IL  , It  and  W 
n 2m  3m  m 

•V? 

W « ^ q (T)  W 

/ . m 

m=  1 


Cl) 


wnere 


+ r;2  q = G (T) 

dT'^  m tn  m 


See  [11].  given  by 


where 


G (T)  = -Ivp.c^a^  f F>'  X,  dX^/N 
m 1 1 ra  1 J n 


2t!p  h a‘*  / (U.?  + + W2)X,  dX- 

s 2m  3m  m 1 3 


p = shell  density,  h = shell  thickness,  P„  = P + P_  + P where  P is  any 

pressure  acting  on  the  shell  which  is  not  applied  through  the  fluid.  For 

a given  shell  geometry  the  in  vacuo  scaled  normal  modes  and  scaled  natural 

frequencies  fi  = aw  /c  may  be  found  by  a variety  of  means,  [11],  and  are 
m m 

therefore  assumed  to  be  known.  The  integral  of  the  unknown  surface  acoustic 
pressure,  P,  in  (23)  may  be  approximated  by 


; P W X.  dX-  - ) H,  P(k)  k'  (k)  X-(k)  (25) 

ml  3 ‘ — k ml 

k 

where  P(k)  = P.  , , i corresponds  to  the  row  of  nodes  for  which  y = 1, 

1 , K. 

k corresponds  to  the  column  number  of  nodes  in  the  y,  plane  and  are  the 
corresponding  weights  in  any  suitable  numerical  quadrature  formula.  Taking 
the  Laplace  transform  of  (22),  assuming  zero  initial  conditions,  gives 


q (s)  = a^G  (s)/(c^(s^  + Q^)) 
m m m 


Equations  (lo),(ll),  (14),  (21),  (25)^  (26)  and  boundary  conditions  are 

suffiritTit  t f>  determine  the  transform  of  the  acoustic  pressure  P , over  the 

1 

region  txtern.il  to  the  shell  and  on  the  sliell  fluid  intcrfate  togetlur  with 
the  transfoni  of  the  normal  shell  deflection.  Solutions  in  the  time  domain 
are  found  thiough  nuitierical  inversion  of  the  Laplace  transform.  See  [12]  and 


[13]. 


Nunerical  Studies 

In  order  to  demonstrate  the  computational  utility  of  the  foregoing, 
numerical  studies  were  carried  out  utilizing  a shell  of  revolution,  i)je 
upper  right  hand  quadrant  of  which  is  shown  in  Fig.  2.  The  scaled  dimensions 
of  the  shell  are  indicated  on  the  figure  together  with  crosses  denoting 
the  approximation  realized  by  retaining  five  terms  in  the  mapping  function 
given  by  (3).  The  mid-section  of  the  shell,  -1  ^ _<  + 1,  is  assumed  to 

be  an  elastic  cylindrical  shell  with  an  unsealed  radius  of  53. 3A  cm.  (21  in), 
thickness  of  2.54  cm  (1  in),  density  of  7.829  • 10^  kg/m’  (G.738  • 10  lb  • 
sec^/in**)  with  a length  of  106.68  cm  (42  in)  for  which  E = 2.0684  • 10^^  Pa 
(30  • 10®  psi)  and  v = 0.3.  The  cylindrical  midsection  is  assumed  to  be 
pinned  at  the  ends  = + 1,  = + 53.34  cm  (+  21  in).  The  end  sections 

of  the  shell,  ] > 1,  are  taken  to  be  fixed  rigid  hemispherical  caps.  The 
shell  is  surrounded  by  an  acoustic  fluid  for  which  c = 1.4478  • 10  m/s 
(0.575  • 10^  in/sec)  with  a density  of  = 0.999  • 10"’  kg/m^  (0.935  • 10 
Ib/sec'^/in  ).  A suddenly  applied  scaled  force,  = cos  (rS^/2)  , acts  on 
the  Interior  of  the  cylindrical  midsection  of  the  shell.  Figure  3 shows  the 
scaled  normal  displacement,  as  a function  of  scaled  time,  at  the  point  = 0 
for  the  shell  surrounded  by  an  acoustic  medium  and  in  vacuo.  Figures  4 and  5 
give  results  similar  to  those  given  in  Figure  3 with  X^  = 0.36726  and  0. 73451 
respectively.  Due  to  the  symmetry  of  the  loading,  only  normal  modes  whirl', 
are  even  functions  of  X^  and  the  corresponding  occur  in  (21)  and  (24). 

Of  these,  it  was  found  necessary  to  sum  only  the  first  two  in  (21)  and  (24). 
The  finite  difference  grid  over  the  first  quadrant  of  the  mapped  aci'usti.- 
field  consisted  of  15  equally  spaced  rows  and  8 unequally  spaced  rclumn.s 
of  nodal  points.  The  nodal  column  spacing  in  the  mapped  acoustic  fiolu  ir 


such  that  the  corresponding  spacing  in  the  physical  plane  is  uniform. 

The  resulting  system  of  equations  was  solved  In  double  precision  for  the 
transforms  of  the  shell  displacements.  Numerical  inversions  were  car-ied 
out  using  [12]. 

Conclus ions 

A numerical  solution  for  the  transient  vibration  of  an  arbitrary  shell 
of  revolution,  surrounded  by  an  acoustic  medium,  has  been  formulated  in 
terms  of  the  conformally  mapped  acoustic  equations,  fluid  shell  boundary 
condition  and  spectral  representation  of  the  normal  shell  displacement. 
Numerical  mapping  techniques  have  been  shown  to  be  effective  in  this  appli- 
cation. The  t ue  readily  includes  the  boundary  condition  at  infinity 

and  is  exac  n those  approximations  implicit  in  the  finite  difference 

method.  1 ired  in  vacuo  normal  modes  and  natural  frequencies  may 

be  computed  utilizing  anyone  of  the  many  dynamic  shell  analysis  codes  now 
available.  Thus  shell  fluid  interaction  problems  for  shells  with  complex 
internal  structure  may  be  treated. 

The  present  formulation,  in  terms  of  the  Laplace  transform,  could  be 
readily  mpdified  so  that  the  time  variable  is  suppressed  through  a finite 
difference  technique.  Steady  state  solutions  may  also  be  derived  through  a 
modification  of  the  boundary  condition  at  infinity. 
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